Theory of magnetic order in the three-dimensional spatially anisotropic 

Heisenberg model 

L. SiurakshinaE and D. Ihle 

Institut fiir Theoretische Physik, Universitdt Leipzig, D-04109 Leipzig, Germany 

R. Hayn 

Institut fiir Theoretische Physik, Technische Universitdt Dresden, D-01062 Dresden, Germany 

(February 1, 2008) 

A spin-rotation-invariant Green's-function theory of long-range and short-range order 
(SRO) in the S = 1/2 antiferromagnetic Heisenberg model with spatially anisotropic 
couplings on a simple cubic lattice is presented. The staggered magnetization, the two- 
spin correlation functions, the correlation lengths, and the static spin susceptibility are 
calculated self-consistently over the whole temperature region, where the effects of spa- 
tial anisotropy are explored. As compared with previous spin-wave approaches, the Neel 
temperature is reduced by the improved description of SRO. The maximum in the tem- 
perature dependence of the uniform static susceptibility is shifted with anisotropy and 
is ascribed to the decrease of SRO with increasing temperature. Comparing the theory 
with experimental data for the magnetization and correlation length of La2Cu04, a good 
agreement in the temperature dependences is obtained. 
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I. INTRODUCTION 

The magnetic properties of spatially anisotropic antiferromagnetic (AFM) quantum spin systems, 
such as the quasi- two-dimensional (2D) parent compounds of high-Tc superconductors (e. g. La2Cu04, 
Ca(Sr)CuO2)0 and the quasi-lD cuprates Sr2Cu03, Ca2Cu03Bi and SrCu02,l are of current interest. 
The main problem is the influence of spatial anisotropy on the staggered magnetization m and the Neel 
temperature T/y in the 3D spin-i AFM Heisenberg model 



(ij)^ {i])y {ij)z 



(1) 



Here Ry — Jy/Jx, Rz = Jz/Jx (throughout we set Jx — 1) and {ij)x.,y,z denote nearest-neighbor (NN) 
bonds along the x-, y- or z-directions of a simple cubic lattice. For real systems, we consider < Rz ^ 

Ry<l. 

In the paramagnetic phase, there exists a pronounced AFM short-range order (SRO) which is re- 
flected by a maximum in the temperature dependence of the magnetic susceptibility at Tmax, where 
0.64 < Tmax < 1.2.i However, the RPA spin- wave theoricsBi and the mean-field theories using auxiliary- 
field representations (Schwinger-bosonJ§j§ Holstein- Primakoff,00 Dyson- Maleev ,00 and boson-fermion 



1 



representationsEJ) which were developed for the quasi-2D model with Ry ^ 1, are valid only at sufficiently 
low temperatures. In those theories, the temperature dependent SRO is not adequately taken into ac- 
count; in particular, the maximum in the magnetic susceptibility cannot be reproduced. In the chain 
mean- field approaches Jil recently improved by spin- fluctuation correctionsEl which lower m and Tjv, an 
asymmetry between intrachain and interchain correlations is introduced. As was shown in Rcf. |^ on the 
basis of a detailed estimate of the exchange integrals for the quasi-lD cuprates using a first-principle 
calculation (Sr2Cu03: Ry ~ 0.004; Ca2Cu03: Ry ~ 0.02), all previous approaches overestimate both 
m and Tn- This deficiency is calling for a theory that provides an improved description of SRO over 
the whole temperature region. In Ref. a spin-rotation-invariant Green's-function theory for the 2D 
isotropic Heisenberg and t-J models was developed, which yields a good description of spin correlation 
functions of arbitrary range and at arbitrary temperatures. Moreover, the susceptibility maximum was 
obtained in good agreement with quantum Monte Carlo calculations. Applying this approach to the 
2D anisotropic Heisenberg model,0 the short-ranged spin correlations at T = are well reproduced as 
compared with exact diagonalization (ED) data. Accordingly, we expect such a theory to describe the 
SRO properties quite well also in the 3D model (|l|). 

In this paper we extend the Green's-function approach of Refs. |l^ and |l^ and present a theory of AFM 
long-range order (LRO) and SRO for the 3D anisotropic Heisenberg model (j^) (Sec. II). Thereby, the 
correlations along all spatial directions are described on the same footing. In Sec. Ill the ground state 
is investigated, where the magnetization and short-ranged spin correlation functions are calculated. In 
Sec. IV we present our finite-temperature results on the i?z-dependence of Tjv, m{T), and of the AFM 
correlation lengths. Moreover, for the first time, the effects of an arbitrary spatial anisotropy on the 
temperature dependence of the uniform static spin susceptibility, especially on T„iax, are investigated. 
The results are compared with experiments on La2Cu04 (magnetization, correlation length, magnetic 
susceptibility) . The summary of our work can be found in Sec. V. 

II. DYNAMIC SPIN SUSCEPTIBILITY 

To determine the dynamic spin susceptibility (Qj'^) = ^(('S'q ! '5'-q))a) by the projection method 
outlined in Ref. we choose the two-operator basis A = {s^, *'^q ) a,nd consider the two-time retarded 
matrix Green's function in a generalized mean-field approximation, ((A;A+))^ — [a; — M'M^^] ^ M 
with M = ([A, A+]) and M' = ([iA, A+]), using Zubarev's notation.^ We get 
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(1) 



(2) 



The spectral moment M^^ = ( 



is given by 



M^'^ = -4Ci, 0,0(1 - cos 9a;) - 4i?yCo,i,o(l - cos^y) - 4^^6*0,0, i(l - cosg^) 



(3) 



The two-spin correlation functions Cr = {Sq ) = Cn^m,i with r = ne^ + mey + le^ are calculated from 

r(i) 
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where n(wq) = (e'^''/'^ - l) . The NN correlation functions are directly related to the internal energy 
per site by e = | (Ci,o,o + RyCo,i.o + RzCo.o,i)- 

To obtain the spectrum in the approximation —S^ = o^qSq , wc take the site representation and 
decouple the products of three spin operators in —S^ along NN sequences introducing vertex parameters 
in the spirit of the scheme proposed by Shimahara and TakadaEi and extending the decoupling given in 
Ref. 0, 

S+S+S[- = a^,'y^^S+Si)S+ + a2{S+Sns; . (5) 

Here af, a\, and af are attached to NN correlation functions along the x-, y- and z-directions, respec- 
tively, and a2 is associated with the longer ranged correlation functions. We obtain 

Wq = 1 + i?^ + i?3 - cosQx - RI cos qy - RI cos 92 + 2alCififi cos{2q^^) - 2a2C'2,o,o cosq^ 
+2Rl{a\CQ^ifi cos{2qy) - a2Co,2,o cos qy) + 2Rl{alCofl,i cos{2q^) - a2Co,o,2 cos g^) 
-2afCififi cosqa; + 2026*2,0,0 - 2i?^(af Co,i,o cos - 026*0,2,0) - 2i?2(aiCo,o,i cosq^ - 026*0,0,2) 
-|-4i?y((af 61,0,0 + "160, 1,0) cosg^; cosg^^ - a26i,i,o(cos g^; -I- cosg^^)) 
-|-4i?2((ai6i,o,o + ai6o,o,i) cosq^ cosq^ - Q!26i,o,i(cosq2; -|- cosq^)) 
+4:RyR^{{a\Cos,() + af 60,0,1) cos cos - £1260,1, 1 (cos +cosg2)) 

-4i?y(ai 61,0,0 cosgj^ + 60, 1,0 cosg^r - 20261,1,0) - 4i?2(ai 61,0,0 cosg^ -I- ai6o,o,i cosg^; - 20261,0,1) 
-4i?yi?2(Q!^6o,i,ocosg2 -f- 60,0,1 cos g^^ - 20260,1,1) . (6) 

Note that our scheme preserves the rotational symmetry in spin space, i.e. x^^i^.^^) = x(Q:'-^) — 
(Qj^)- For |q| ^ 1 we have 



„2 „2 , 2 „2 



„2„2 



= c^Qx + Cyqy + c^g^ 



(7) 



with the squared spin-wave velocities 

= 2 ~ 3aiCi,o,o + 0^2^2,0, - 2Ry{aiCifi,o - a2Ci,o,o) - 2Rz{aiCifi,o - a2Ci,o,i) , (8) 

= Rli^ - 3a5'Co,i,o + a2Co,2,o) - 2Ry{a\Co,i.o ~ a2Ci,i,o) - 2RyR^{a\CQ.ifi - q;2Co,i,i) , (9) 

and 

cl = Rl{\ - ialCo^o,! + a2Co,o,2) - 2i?,(a^Co,o,i - ^2^1,04) - 2i?,yi?,(a^Co,o,i - ^2^0,14) • (10) 

Considering the uniform static spin susceptibihty x = l™q^o -'^^q^V(2'-^q)i the ratio of the anisotropic 
functions and Wq must be isotropic in the hmit q ^ 0. That is, the conditions 

[cy/c^f = RyC^Afl/Ci^ofi (11) 

and 

(c,/c,)2 = i?,Co,o,l/C^l,0,0 (12) 

have to be fulfilled. 

The critical behavior of the model (|^) is reflected in our theory by the closure of the spectrum gap at 
Q = (tt, tt, tt) as T approaches T^r from above, so that limr^Tiv X~^(Q) — 0- At T < Tj^ we have ujq — 
and, separating the condensation part C, 

Cr = ^Y. Cqe''"" + Ce'Q^ (13) 
q(#Q) 

where C results from ( p^ with r = employing the sum rule Co. 0,0 = \- Then the staggered magneti- 
zation m is calculated as 

™' = ^E(SoSr)e-^"- = ^C. (14) 

r 

The theory has 14 quantities to be determined self-consistently (9 correlation functions in w^, m, and 
4 vertex parameters) and 13 self-consistency equations (10 Eqs. ( |T^ ) including Cq.o.o = \: the LRO 
condition ljq ~ 0, and Eqs. ( p4| ) and (p^)). If there is no LRO, we have ljq > 0, and the number 
of quantities and equations is reduced by one. As an additional condition for determining the free a 
parameter at T = 0, we adjust the ground-state energy per site which we compose approximately as 
e(i?y, Rz) — t{Ry; 0) + e(0, Rz) ~ e(0, 0), where e{Ry, 0) (and e(0, Rz)) is taken from the Ising-expansion 
results by Affleck et al. for the 2D spatially anisotropic Heisenberg model,0 and e(0, 0) = —0.4431 is the 
Bethe-ansatz value. This approximation is suggested to be good at least for Rz ^ Ry (or Ry -C Rz)- 
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To get an additional condition also at finite temperatures, where e data are not available and all vertex 
parameters are temperature dependent, we assume, following Refs. and p^, the ratio 

as temperature independent. 

III. GROUND-STATE PROPERTIES 

In Fig. 1 our results for the zero-temperature staggered magnetization ttlq = m{T = 0) as a function 
of Ry and Rz are shown. They indicate an order-disorder transition at the phase boundary Rz.c{Ry) 
or Ry.c{Rz) (cf. inset). For i?^ = we get the critical ratio Ry^dO) — 0.24 which was already found 
m Ref. In that paper the suppression of LRO below the finite value of Ry^c was interpreted, in 

combination with ED data, as indication of a rather sharp crossover in the spatial dependence of the spin 
correlation functions in the LRO phase at the coupling ratio Ry_o ~ 0.2. The finite value of Ry^c, however, 
seems to be due to the approximations in our theory, since there are strong indications for Ry,c = (see 
Ref. Accordingly, we cannot explain the tiny magnetic moments of Sr2Cu03 and Ca2Cu03j^ since 
for Ry ^ l,i we have m = 0. This result is just opposite to the overestimation of m by all previous 
spin- wave theories. As seen in the phase diagram (inset of Fig. 1), the inclusion of the interplane coupling 
Rz stabilizes the LRO, where this effect is quite considerable even at very small values of Rz- 

Figure 2 exhibits some short-ranged spin correlation functions at T = 0. For Rz = 0, in Ref. |l^ the 
correlators Ci,o.Oj Co.i,Oi a-nd Ci.i^o a-s functions of Ry were found to agree well with the ED data. For 
Rz = 0.02 (cf. Fig. 2) our results deviate only slightly from those at Rz = 0. The sign changes and 
magnitudes of Cr reflect the AFM SRO. In the limit Ry the correlations between the x-z-planes 
vanish. At Rz > Rz.c{0) — 0.24 the LRO enhances the inter-x-z plane correlators and results in their 
sharp drop towards their limiting value Ce'*^"" as Ry 0. This is visible in the data for Rz = 0.35 in 
Fig. 2. 

IV. FINITE-TEMPERATURE RESULTS 

At nonzero temperatures we have solved the self-consistency equations (|l^) supplemented by the con- 
ditions (pT[), (p^, and (|l5|) to obtain the magnetization m{T), the Neel temperature {m{TN) — 0), the 
static spin susceptibility, and the anisotropic correlation lengths. 
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In Fig. 3 the Neel temperature is plotted as a function of . For i?^ = we get T^v = (see Ref . ^6|) , 
in agreement with the Mermin- Wagner theorem. The increase of T/y with Rz is governed by the intra-x-y 
plane anisotropy. At a fixed value of R^ , the decrease of T/v with decreasing Ry is in accordance with the 
reduced zero-temperature magnetization (cf. Fig. 1). Comparing our results for Ry — 1 with previous 
RPA/mean- field approaches (see Table I), we ascribe the reduction of T/v as compared with Refs. ^||, 
and ^ to the improved description of SRO. That is, the LRO is suppressed in favor of a paramagnetic 
phase with pronounced AFM SRO. If Rz is fit to the Neel temperatures of real systems, the strong 
overestimation of T/v by previous theories results in very small values of the interplane coupling. In our 
approach the resulting Rz values turn out to be higher. Considering La2Cu04 with T/v — 325Ki0 and 
putting J = ISOmeV {J = = Jy) or J = 117meV,lli we obtain Rz ~ 10"''' or Rz ~ 1.6 x 10-^ 
respectively, in contrast to Rz < 10""* according to Refs. || and For Cao.85Sro.i5Cu02 (T/v = 540K, 
J = 125meV)il we get i?2 ~ 1.2 X 10 ^ as compared with Rz ~ 2.5 x 10 ^ obtained from a fit of the 
low-temperature magnetization data.@ 

Figure 4 shows the temperature dependence of the staggered magnetization at Ry = 1 (for the zero- 
temperature values, compare with Fig. 1). The shape of the normalized curve m/mo versus T/Tn (see 
inset) depends on the single parameter Rz and is similar to that found in previous spin-wave theories.B~E0 
At low enough temperatures the system exhibits 3D behavior, so that the decrease of m follows a T^ law. 
This was also observed by NMR experiments on LasCuoJi] (T/v = 312K) yielding m/mo = l-a{T/TNf 
with a — 0.67 for T < lOOK. The NMR data is indicated in the inset of Fig. 4 (marked by a bold curve) 
and agrees well with our theory for Rz = 10"'' (as estimated above). For temperatures close to T/v our 
numerical results for m{T) are described by the law m{T) cx (1 — T/TnY^^. The square-root temperature 
behavior agrees with the findings of Refs. [7|j8|JlO|, and with the neutron scattering data on La2Cu04,i but 
contradicts the result of Ref. |l^ (m oc 1 — T/T/v). 

Considering the AFM correlation lengths above T/v and for Ry ^ 1, the expansion of x(q) around Q, 
x(q) = x(Q) [l + ^xyi^x + ^y) + Cz^z] ^ ^^th k = q — Q, yields the intraplane correlation length 

'1 



-2 



Q 



2 + ll<Ci^o,o + «2 (Ca^o.o + 2Ci4,o) + 



2C 

+2Rz (a^Ci.o.o + 2a^Co,o,i + «2Ci,o,i)] 7^ (16) 



and the interplane correlation length 

£f = -RzOJ^^ [4 (2a^Ci,o,o + a?Co,o,i + a2Ci,o,i) + 
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. (17) 



In Fig. 5 the influence of the interplane couphng on the temperature dependence of S^^y and ^ (inset) 
is shown. For comparison, the intraplane correlation length at i?^ = (see also Ref. ^6|) is plotted, where 
the low-temperature expansion £,^y = 2(2af |Ci^o,o(0)|)^/^T-i exp[27raf r7i^/(3r)] holds up to T = 0.2 
within a deviation of about 6% from the full temperature dependence calculated by Eq.(|l^). For i?^ > 
the correlation lengths diverge at Tjv, since the gap wq closes as T approaches Tjv from above. In the 
vicinity of Tjv, S,^y and behave as T — Tjv also found by previous mean-field approaches El'El 

Let us compare our results for the intraplane correlation length with the neutron-scattering data on 
La2CuO40 in the range 340K< T < 820K shown in Fig. 6. Taking J as obtained previouslyEl from a 
least-squares fit of ^^y in the 2D model (a = 3. 79 A), J = 117mcV, for T > 500K and i?^ < 3.5 x 10"^ 
get a good quantitative agreement with experiments. In Ref. |l^ the deviation of the theory for Rz — 
and T < 500K from the experimental data was ascribed to the appearance of the preexponential factor 
T^^ in the low-temperature expansion of £^xy which is an artifact of our mean-field approach. However, 
this deviation may be reduced by the inclusion of the interplane coupling, since ^^yC^N) = 0. For 
Tjv — 325X0 we get ~ 1.6 x 10"'^ (see above. Fig. 3), and the theoretical ^~y curve lies between 
the Rz — result and the experiments. The discrepancy between the theoretical and experimental low- 
temperature correlation lengths may be further reduced by the choice of higher Rz values. Taking, for 
example, Rz = 3.4 x 10^"^, we get a very good quantitative agreement (cf. Fig. 6) down to 360K; however, 
the Neel temperature turns out to be somewhat higher (T^v — 353K). 

Finally, we consider the uniform static spin susceptibility x(T) = limq^o x(ci)- In Fig- 7 the anisotropy 
effects on the temperature dependence are demonstrated. For Rz — and a strong intraplane anisotropy 
{Ry < 0.2) the minimum of x{T) at a finite temperature, being an artifact of our approach, may signal 
the crossover in the spatial dependence of the spin correlation functions at i?y,o — 0.2 as was discussed in 
Sec. III. Note that such a minimum in the ID model {Ry = 0) was also found in Ref. At Ry > 0.2, 
the increase of x with temperature, the maximum at Tmax near the exchange energy = I (see inset), 
and the crossover to the high-temperature Curie- Weiss behavior are due to the decrease of AFM SRO 
with increasing temperature (cf. Ref. p^ ). Let us point out that the susceptibility maximum is totally 
missed in RPA theories.0 With increasing Ry, we obtain an increase of T„iax which agrees with a general 
tendency found in various spin-1/2 Heisenberg models and analyzed in Ref. |^. For comparison, the exact 
values at Ry = and Ry = 1 are given by T^ax ~ 0.640 and T^ax = 0.94,0 respectively. Since our 
theory allows the calculation of Tmax at any spatial anisotropy, it may provide a reliable interpretation 



of experimental data on low-dimensional spin systems. Considering the maximum spin susceptibility 
Xmax — x{Tmax), again our results are in accordance with the general behavior:@ Xmax increases with 
decreasing T^axj i-S- with decreasing Ry. Concerning the influence of the interplane coupling, the 
enhancement of the low-temperature susceptibility by Rz may be explained by the weakening of the SRO 
effect in higher dimensions. As seen from Figs. 7 and 3, the uniform susceptibility has no singularity at 
the Neel temperature, contrary to the RPA result of Ref. revealing a peak of x{T) at Tjv. Concerning 
the maximum in xiT) of La2Cu04, we get T^ax = 1.19 J = 1615K (cf. Fig. 7, J = 117meV). This value 
roughly agrees with the estimate given by JohnstonJ^^ Tmax = 1460iir, by means of a scaling analysis of 
the susceptibility data below 800K. 

V. SUMMARY 

In this paper we have extended the spin-rotation-invariant Green's-function theory of magnetic LRO 
and SRO in 2D Heisenberg modelJlllli to the 3D Heisenberg antiferromagnet with arbitrary spatial 
anisotropy. Our theory provides a satisfactory interpolation between the low-temperature and high- 
temperature behavior, where the temperature dependent SRO, described in term of two-spin correlation 
functions, is adequately taken into account. The main results are summarized as follows. 

(i) The incorporation of SRO results in a strong suppression of Neel order with increasing anisotropy 
and in a reduced Neel temperature as compared with previous spin-wave approaches. 

(ii) The temperature dependence of the uniform static spin susceptibility reveals a maximum in the 
short-range ordered paramagnetic phase and a crossover to the Curie- Weiss law. The position of the 
maximum is influenced by the spatial anisotropy. 

(iii) Comparing the theory with experiments on the magnetization and correlation length of La2Cu04, 
a good quantitative agreement is found. 

From the results of our theory we conclude that the application of this approach to extended Heisen- 
berg models (anisotropy in spin space, frustration) may be promising to describe the SRO effects on the 
unconventional magnetic properties of real low-dimensional spin systems. 
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TABLE I. Neel temperature T^/ aX Ry = \ compared with other approaches 
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Figures 

Fig. 1. Staggered magnetization at T = as a function of spatial anisotropy. The inset shows the 
stabihty region of Neel order. 

Fig. 2. Spin correlation functions at T = for different spatial anisotropies. 
Fig. 3. Neel temperature as a function of ~ Jz/ Jx- 

Fig. 4. Staggered magnetization vs. temperature for Ry = 1. The inset shows the R^ dependence of 
the normalized curves compared with the NMR data on La2CuO40 (bold curve). 

Fig. 5. Inverse antiferromagnetic correlation lengths within {i,^y) and between the x-y planes (■fj^, see 
inset) for Ry — 1. 

Fig. 6. Inverse antiferromagnetic intraplane correlation length in La2Cu04 obtained by the neutron- 
scattering experiments of Ref. ^ and from the theory {Ry — 1) for different Rz values. 

Fig. 7. Uniform static spin susceptibility vs. T. The inset exhibits the position Tmax of the maximum 
in x{T) vs. Ry. 
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